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ABSTRACT 

We revisit the Blandford & Znajek (1977) process and solve the fundamental equation that governs 
the structure of the steady-state force-free magnetosphere around a Kerr black hole. The solution 
depends on the distributions of the magnetic field angular velocity w and the poloidal electric current 
/. These are not arbitrary. They are determined self-consistently by requiring that magnetic field lines 
cross smoothly the two singular surfaces of the problem, the inner 'light surface' located inside the 
ergosphere, and the outer 'light surface' which is the generalization of the pulsar light cylinder. We 
find the solution for the simplest possible magnetic field configuration, the split monopole, through a 
numerical iterative relaxation method analogous to the one that yields the structure of the steady-state 
axisymmetric force- free pulsar magnetosphere (Contopoulos, Kazanas & Fendt 1999). We obtain the 
rate of electromagnetic extraction of energy and confirm the results of Blandford and Znajek and of 
previous time dependent simulations. Furthermore, we discuss the physical applicability of magnetic 
field configurations that do not cross both 'light surfaces'. 

Subject headings: Physical Data and Processes: Accretion, accretion disks — Black hole physics — 
Magnetohydrodynamics (MHD) 
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1. INTRODUCTION 

The electromagnetic extraction of energy from a spin- 
ning black hole can be one of the most powerful and ef- 
ficient engines in the Universe. It was Blandford & Zna- 
jek (1977) (hereafter BZ77) who first argued that a spin- 
ning black hole threaded by a sufficiently strong magnetic 
field and permeated by an electron-positron plasma gen- 
erated from pair cascades, establishes a force-free mag- 
netosphere, in direct analogy to the theory of the ax- 
isymmetric pulsar magnetosphere developed a few years 
earlier by Goldreich & Julian (1969). They derived the 
fundamental equation that governs the structure of this 
magnetosphere, the general relativistic force-free Grad- 
Shafranov equation, and obtained a second order pertur- 
bative solution for small values of the Kerr parameter 
a and a scaling for the electromagnetic power extracted 
from the rotating black hole. 

The structure, singular surfaces and general prop- 
erties of this equation for magnetic fields threading 
Schwarzschild and Kerr black holes were discussed re- 
cently by Uzdensky (2004, 2005). He emphasized the 
need for a self-consistent determination of the distri- 
butions of the magnetic field angular velocity and the 
poloidal electric current, something that early relaxation 
methods (Macdonald 1984) and later higher order per- 
turbative solutions (Tanabe & Nagataki 2008) did not 
take into account. He also noted that this may be a 
very difficult task, and indeed, to the best of our knowl- 
edge a solution of this equation in the general case of 
an open magnetic field configuration threading a Kerr 



black hole is still lacking. Nevertheless, the broader 
astrophysical interest in the Blandford-Znajek process, 
has prompted in the past decade several researchers to 
seek and obtain the steady-state structure of the spin- 
ning black hole magnetospheres through time-dependent 
general relativistic force-free numerical simulations (e.g. 
Komissarov & McKinney 2007, Palenzuela et al. 2011, 
Lyutikov & McKinney 2011). These works established 
the universal applicability of the electromagnetic extrac- 
tion of the rotation energy of a Kerr black hole in pow- 
ering active galactic nuclei, jets in X-ray binaries, and 
even gamma-ray bursts. 

The presence of this extra source of energy in addition 
to that of the accretion disk may account for some of 
the phenomenology. One should bear in mind, however, 
that the situation is more complicated. For instance: (1) 
Almost 90% of extragalactic sources are radio quiet, i.e. 
they do not produce electromagnetic outflows in the form 
of radio jets, and yet they are believed to contain both 
constituents of the Blandford-Znajek process, namely a 
spinning black hole and a strong magnetic field thread- 
ing it (e.g. Kukula 2003); (2) The formation and disrup- 
tion of radio jets in black hole X-ray binaries over short 
time scales is at odds with this general notion (e.g. Bel- 
loni 2010); (3) There seems to exist strong observational 
evidence that the power in the radio jet is not directly 
related to black hole spin, contrary to the main result of 
BZ77 (Fender, Gallo & RusseU 2010; however, see also 
Narayan & McClintock 2012). Notice, though, that this 
may be due to an erroneous determination of the black 
hole spin. 
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The above justify an independent re-evaluation of the 
Blanford-Znajek process through the solution of the gen- 
eral relativistic force-free Grad-Shafranov equation. The 
special relativistic form of that equation, the so called 
pulsar equation, was only solved in 1999 by Contopoulos, 
Kazanas and Fendt who showed that the magnetic field 
structure can only be obtained together with the deter- 
mination of the unique electric current distribution that 
allows for a smooth and continuous structure to exist. 
In that respect, steady-state solutions offer a deeper in- 
sight in the physics of the problem than numerical time- 
dependent ones. It should be noted that, while BZ77 
provided an estimate of the energy extraction efficiency 
of this process, they did not obtain an exact solution of 
the structure of the black hole magnetosphere. Also, at 
that time, the critical role of the singular surfaces and 
the distributions of the magnetospheric electric current 
and field line angular velocity had not been appreciated. 
Today, we revisit this problem with all the knowledge we 
carry from our 13-year long investigation of the force-free 
pulsar magnetosphere. 

2. THE GENERAL RELATIVISTIC PULSAR EQUATION 

In order to derive the fundamental equation that gov- 
erns the steady-state structure of the force- free magneto- 
sphere around a Kerr black hole we follow closely the 3-1-1 
formulation of Thorne & Macdonald (1982) used by most 
researchers in the astrophysical community (e.g. Uzden- 
sky 2005). We restrict our analysis to steady-state and 
axisymmetric space times where (. . .)^t = (■ • ■),<j> = 0. In 
that case, the general 4-dimensional space-time geometry 
may be written in Boyer-Lindquist spherical coordinates 
x'^ = {t, r, 9, (j)) as 



ds^ = g^iidx^dx'^ 



-a^dt^ 
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Here, a = (AS/A)!/^ ^nd fl = 2aMr/A are the lapse 
function and angular velocity of 'zero- angular momen- 
tum' Hdncial observers (ZAMOs) respectively. 



A 



= r2 



2Mr 



S = + cos^ 9 



A= (r^ +a^)^ -a^Asin^6l , 



M is the black hole mass, and a its angular momentum 
(0 < a < AI). Throughout this paper we adopt geo- 
metric units where G = c = 1. Semicolon stands for 
covariant derivative, comma for partial derivative. Latin 
indices denote spatial components (1 — 3), Greek indices 
denote space-time components (0 — 3), and '~' denotes 
the spatial part of vectors. ZAMOs move with 4- velocity 
= (1/a ,0 ,0 ,ri/a) orthogonal to hypersurfaces of 
constant t. The force-free magnetosphere of a spinning 
black hole is characterized by the electromagnetic energy- 
momentum tensor 



(2) 



and the condition T.^"^ — 0. Here, the rest mass 
and pressure contribution have been neglected. The 



electromagnetic field tensor F'^'^ is related to the elec- 
tric and magnetic fields , i?^ measured by ZAMOs 
through Ft'" = VE" - VEt" + ef'^^PExUp {e^^xp = 
[/ii^Ap]|det(g^y)|^^/^ is the 4-dimensional Levi-Civita 
tensor). Under these conditions, the fundamental equa- 
tion that governs the steady-state structure of the force- 
free magnetosphere around a Kerr black hole becomes 



PeE + JxB^O. 



(3) 



Pe and J are the electric charge and current densities re- 
spectively. Eq. ([3|) is supplemented by Maxwell's equa- 
tions of electrodynamics 

V -B^O 

V ■ E = 47rpe 

V X (aB) = 4:Tra J 

V X («£■) = . 



(4) 



Here, 



V • = Af^. , A-B = g'^A.B, , 
(Vxi)^ = [zjfc]|det(.g„„) 1-1/2^,;, , 

(i X By = [ijk]\det{gim)\-^'^A,Bu 

For several applications in astrophysics, perfect (infinite) 
conductivity is a valid approximation. In this case, 

E-B = 0, (5) 

and the electric and magnetic vector fields can be ex- 
pressed in terms of three scalar functions, ^{r, 9), uj{'^), 
and as 



B{r,9) 



1 



y/Asm9 
fl — UJ 



*,e,-VA*,„^— ^ (6) 



E{r,9)^^^{VA^,r,^,0,o} 



(7) 



u is the angular velocity of the magnetic field lines, / is 
the poloidal electric current flowing through the circular 
loop r =const., 9 =const., and VE* is equal to (Stt)"^ times 
the total magnetic flux enclosed in that loop. Notice that 
the electric field changes sign close to the horizon with 
respect to its sign at large distances. As explained in 
BZ77, a rotating observer (ZAMO) will in general see 
a Foynting flux of energy entering the horizon, but he 
will also see a sufficiently strong flux of angular momen- 
tum leaving the horizon. That ensures that energy is 
extracted from the black hole. The poloidal component 
of Eq. ([3]) then yields the general relativistic force-free 
Grad-Shafranov equation 
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(Eq. (3.14) of BZ77 re- written in our notation). Hence- 
forth, primes will denote differentiation with respect to 
5". One sees directly that if we set a = and M = in 
eq. ([2]) we obtain 
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1 cos 6* 



[1- 



cj^r^ sin^ ( 



2*, 



- 2uj^ cos sin 6*^1^ , 



- ujuj'r^ sin^ 9 {^^^^ + j = -4//' , (9) 

which is the well known pulsar equation (Scharlemann & 
Wagoner 1973). The zeroing of the term multiplying the 
second order derivatives in eq. 

Lo'^Asi-a^e AMaujrsm^e 2Mr , , 

1 S ^-0' (^0) 

yields the singular surfaces of the problem. We will 
henceforth call the singular surfaces 'light surfaces' 
(LSj3. When we set M = a = 0, it yields the stan- 
dard pulsar light cylinder rsinO = c/uj. In general, the 
shape of the outer LS is only asymptotically cylindrical 
as 6* — (see figure 2 below), and an inner LS appears 
inside the ergosphere. We will henceforth use the no- 
tation 'light cylinder' (LC) and 'outer LS' interchange- 
ably. It is interesting to note that the outer boundary of 
the ergosphere corresponds to the solution of the singu- 
larity condition for a; = 0, whereas the inner boundary 
(the event horizon at r = tbh = M + y/lVP — a?) cor- 
responds to the solution of the singularity condition for 
u) = 17bh = a/(r|jj + a^), where 17bh is the angular ve- 
locity of the black hole. It is also interesting to note that 
the natural 'radiation condition' at infinity (energy must 
flow outwards along all field lines) requires that 

< w < 17bh (11) 

(BZ77), and therefore indeed the inner LS lies inside the 
ergosphere. 

Both Eqs. ^ and © contain the two functions, cj(^) 
and which must be determined by the physics 

of the problem. In the case of an axisymmetric spin- 
ning neutron star, uj is usually taken to be equal to the 
neutron star angular velocity ^InA In pulsars, /(^) 

^ These are none other than the 'velocity- of-light surfaces' of 
Macdonald (1984) where the speed of a particle moving purely 
toroidally with angular velocity u! equals c. 

^ Notice that in the presence of particle acceleration magne- 
tospheric 'gaps', this is not 100% exact (Ruderman & Suther- 
land 1975, Contopoulos 2005). In particular, in old pulsars near 
their death line uj S^ns- 



is self-consistently determined through an iterative nu- 
merical technique that implements a smooth crossing of 
the relativistic Alfven surface, which is the usual light 
cylinder where rsinO = c/uj (Contopoulos, Kazanas & 
Fendt 1999, Timokhin 2006). In the case of a spinning 
black hole, the situation is qualitatively similar but more 
complicated. Contrary to a neutron star, the black hole 
does not have a solid surface, and therefore it cannot im- 
pose any restriction on lo other than the 'radiation con- 
dition' Eq. ([TT|) . The only natural restriction that is im- 
posed by the physical problem is that magnetic field lines 
must be smooth and continuous everywhere. a;('I') must 
therefore be determined together with /(^) through the 
condition of smooth crossing of both LS of the prob- 
lem, the inner one inside the ergosphere, and the outer 
one at large distances. In the next section, we will see 
that this can be achieved through an iterative numerical 
technique analogous to the one employed in the pulsar 
magnetosphere. As remarked previously, iterating with 
respect to two functions simultaneously is indeed a very 
difficult task. This is why Uzdensky (2005) opted to con- 
sider only field lines that connect the black hole horizon 
to the surrounding accretion disk where a;('I') is deter- 
mined by the Keplerian angular velocity of the disk. 

There exists one more interesting complication. 
Eq. pri) applies only to open field lines that cross the 
event horizon. In the case of a star, field lines that cross 
the stellar surface rotate with the angular velocity of the 
star, but they do not cross an event horizon (either be- 
cause no event horizon forms, or because they simply 
avoid the horizon), and the above restriction does not 
apply. In the present case, we are interested in studying 
the fundamentals of the Blandford-Znajek process, and 
therefore, we will consider a black hole 'as clean as possi- 
ble'. Obviously, a magnetized astrophysical black hole is 
not isolated. In that respect, we would like here to study 
an astrophysical black hole with the minimum number of 
extra elements, namely a surrounding thin disk of matter 
to hold the magnetic field through the necessary electric 
currents and charges. Charges and currents will be pro- 
duced through pair production in the black hole magne- 
tosphere. We will ignore magnetic field lines that do not 
cross both the inner and outer LS because we will have 
no physical way to determine their to and the poloidal 
electric current / contained inside them. We will also 
ignore magnetic field lines that cross the equator out- 
side the black hole event horizon because such field lines 
need a source of poloidal electric current on the equa- 
tor (if they rotate, they will need to cross at least one 
LS where the condition of smooth crossing will in general 
require a nonzero electric current to flow along that line) . 

Therefore, in the present work we will only study black 
hole magnetospheres where all magnetic field lines reach 
the event horizon, and all magnetic field lines are open. 
As argued in Lyutikov (2012), such configuration is the 
most natural end stage of stellar collapse. Notice once 
again that it is sustained by electric currents in an equa- 
torial current sheet. The effect of a surrounding disk of a 
certain height possibly threaded by the return magnetic 
field (as in the scenario predicted by the Cosmic Bat- 
tery; Contopoulos & Kazanas 1998) will be considered 
in a future work. 
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3. THE NUMERICAL ITERATIVE SCHEME 

We will here describe the technical numerical details 
that led to the solution of Eq. ([2]). Firstly, we change 
the radial variable from r to R{r) = r/{r + AI). Our 
numerical integration extends from i?niin = -R(?'bh) (the 
event horizon) to i?max = 1 (radial infinity). The 9 co- 
ordinate extends from 0min = (the axis of symmetry) 
to ^niax = 7r/2 (thc equatorial plane). We implemented 
a 256 X 64 numerical grid uniform in R and 9. Notice 
that this grid has a very high resolution in r around the 
black hole event horizon where the inner LS lies, but not 
as high around the outer LS, and in particular at low 9 
and a values. 

Secondly, we specify boundary conditions on the axis of 
symmetry, the horizon, the equatorial plane, and infinity. 
^'(6' = 0) = 0. ^'(r'BH, ^) is determined by the integration 
over 9 of the so called Znajek regularization condition 



-0.1[*(r+,e') - *(r-,e')] , (15) 



/(*) 



MrBH sin ^ 



-(17bh - w)^',, 



(12) 



' BH 



(Znajek 1977). Notice that this condition is updated as 
the distribution 7(5') is updated along with the numeri- 
cal solution for '^{r,9). The equatorial boundary condi- 
tion is interesting. As we discussed in the previous sec- 
tion, the equatorial region may contain an accretion disk, 
with or without its own magnetic field. Before investigat- 
ing such a complex astrophysical system, we would like to 
understand first the simplest case, that of a magnetized 
black hole where the electric currents supporting its mag- 
netic field are distributed on a thin equatorial current 
sheet. In that case, the equatorial boundary condition 
becomes \l/(r, = 7r/2) = VPmax- The outer boundary 
condition is also not important since the outer boundary 
has practically been moved to infinity. In practice, we 
chose ^,r(^ = 1) = 0, but any other boundary condition 
yields the same results. We initialize our numerical grid 
with a split monopole field configuration 

*(r, 9) = ^suir, 9) = ^',„ax(l - cos 9) , (13) 



with 



a;(*) = 0.5f2BH and 



/(*) - /sm(*) = -0.25f)BH*[2 - {^/-^nrax)] ■ (14) 

The reader can directly check that Eqs. & (|14l) are 
an exact analytical solution of the pulsar equation (Eq.[9l 
Michel 1982). Notice the negative sign in Eq. For 
the particular magnetic field configuration where on 
the axis is aligned with the black hole spin direction, 
this corresponds to an outflow of electrons. Different ini- 
tial configurations yield the same final solution. We dis- 
cretize all physical quantities on our 256 x 64 {R, 9) grid 
and we update '^{R,9) through simultaneous overrelax- 
ation with Chebyshev acceleration (subroutine SOR from 
Numerical Recipes; Press, Flannery & Teukolsky 1986). 

Thirdly, we update the distributions of uj{'ii) and /(^) 
as follows: At each latitude 9, we check where the sin- 
gularity condition fEgfTO)) is satisfied in r. At each such 
radial position, we extrapolate ^E* inwards from larger r 
(5'(r+,6')) and outwards from smaller r (^(r^,6')). In 
general, 'i/{r~^,9) and 5'(r~,0) differ. Then, at the inner 
LS we implement 

^new (^'new) — Wold(^' new^ 



whereas at the outer LS we implement 

-^new(^'new) = -^old (4'iiow) 

+0.1[*(r+,e') - *(r",( 

(^'now) =1^01(1 (* new J 

-0.1[^'(r+,6') - *(r",( 



where 



^'ncw = O.5[^'(r+,0) + *(r-,( 



(16) 
(17) 



at each LS. The reasoning here is that we impose 
weighted corrections on uj{^) and 1(5') based on the 
non-smoothness of the \I'(r, 9) distribution along all grid 
points inside and outside the two LS, and not through 
the regularization conditions as we did in Contopoulos, 
Kazanas & Fendt (1999). This is a very general proce- 
dure that may be applied to any similar singular equa- 
tion. 

Our numerical scheme does not converge easily, and 
its parameters need to be empirically adjusted by fol- 
lowing the progress of the iteration. The signs of the 
weighting coefficients are determined by trial and er- 
ror: when we choose the wrong signs, the mismatches 
\'${r^ ,9) — 'i>{r~ ,9)\ increase and the^'(r,0) distribution 
is very quickly destroyed. Their magnitudes are also ob- 
tained empirically: too large values (around unity) lead 
to over-corrections and instability, too small values (be- 
low 0.01) do not correct fast enough for the iteration 
to converge. Furthermore, in order to facilitate conver- 
gence, at the inner LS we update only w, whereas at 
the outer LS we update both / and lu every other relax- 
ation iteration for 4*. Also, in order to avoid numerical 
instabilities, we smooth out the distributions of / and 
u! every 50 relaxation iterations for v]/. We found that 
more frequent smoothing inhibits convergence. Our goal 
is to minimize the residuals in the discretized form of 
eq. ([2]) together with the non-smoothness in 'i'{r,9), but 
our method fails for poor grid resolutions at the inner 
and/or the outer LS (see Discussion section). 

The results of our numerical integration are shown in 
figure 1. We show here the distributions of aj(^) (top 
panel), 2|/(5')| (center panel), and Foynting fiux at large 
distances 2|/(^)|ti;(5') (bottom panel) after 6000 relax- 
ation steps for '^{R,9), for six values of the black hole 
spin parameter a. As is the case in pulsars, the mag- 
netospheric electric current is non-zero at = ^'max, 
and the global electric circuit closes through an equa- 
torial current sheet. Notice that we have no way to 
update uj along the axis, and therefore we have cho- 
sen uj(JSj = 0) = 0.5r2BHEI On the equator we imple- 
mented uj'{^ = "iimax) = /'(^' = I'ma:^) = 0. The reader 
should check the qualitative agreement with the results 
of the time-dependent force-free general relativistic sim- 
ulations performed by Tchekhovskoy, Narayan & McK- 
inney (2010) (their figure 4). 

In figure 2 we show the characteristic magnetic field 
configuration for a/M = 0.9999 near the black hole hori- 
zon (left panel) and at larger distances (right panel). 
We also plot the inner LS inside the ergosphere, and 

^ We have also run our simulations with u)'{'i! = 0) = and 
the solutions were similar, only the convergence of ui became more 
unstable. 
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Figure 1. Distributions of lj('I') normalized to SIbh (top panel), 
2|/(*)| normalized to SIbh* max (center panel), and Poynting flux 
at large distances 2|7(^)|a;('I') normalized to Hgjj'I'maa: (bottom 
panel) for various values of the black hole spin parameter a/M = 
0.7 (solid), 0.8 (short dashed), 0.9 (dotted), 0.99 (small dots), 0.999 
(dash dotted), 0.9999 (short dash dotted). All three distribut ions 
are close to the analytical split monopole expressions of Eq. H14| l 
(thick dots). 



the outer LS which becomes asymptotically cylindrical 
as — >■ 0. Notice that magnetic field lines are monopo- 
lar at large distances (beyond a few times the radius of 
the event horizon), and they bend upwards towards the 
axis as they approach the horizon. This effect is also 
seen in the numerical simulations of Lyutikov & McKin- 
ney (2011) for a = 0.99. It is not discernable in the field 



Figure 2. Poloidal magnetic field lines (lines of constant 't in 
20 equal subdivisions of '^max) for a/M = 0.9999 close to the 
event horizon (top panel) and further out (bottom panel). The 
thick lines are the event horizon, the inner LS, the outer boundary 
of the ergosphere, and the outer LS (from the center outwards 
respectively). 

configurations of figure 1 of Tchekhovskoy, Narayan & 
McKinney (2010) because the scale is not adequate, but 
is implied in their figure 7. Notice that is equal 

to 0.5i7BH to within 10%, and that /(\I/) is very close to 
the split monopole expression /sm(^) also to within 10% 
(Eqs.[Tn]&[Tl]). Therefore, the Znajek condition (Eq-HT 
can be rewritten using Eq. ^ as 



2MrBH 



(2M)2 



{a/Mf 



[1 + v/1 - HM f 



■ cos^ e 



(18) 



Finally, we calculate the total rate of electromag- 
netic energy extraction (from both hemispheres) as mea- 
sured at infinity through Eq. (4.11) of Blanford & Zna- 
jek (1977) as 



£: = -2 X 27r 



2/(*)w(*)d^' 



27r, 



(19) 



We have here normalized our results to the split 
monopole value with wsm = O.SrisH- k sa 0.8 for 
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a/M = 0.9999 and is practically indistinguishable from 
unity for a/M <, 0.9 (see bottom panel of figure 1), thus 
confirming the main result of BZ77, namely that, in an 
open magnetic field configuration with an infinitely thin 
surrounding equatorial disk, electromagnetic energy is 
extracted at a rate proportional to both ri|jj and '^^ax- 

4. DISCUSSION 

We believe that the study of the steady-state open 
force-free black hole magnetosphere offers a deeper in- 
sight in the physics of the problem than time-dependent 
numerical simulations. In analogy to the axisymmetric 
force-free pulsar magnetosphere, we argued that steady- 
state solutions can only be obtained through a self- 
consistent determination of the distributions of both the 
magnetic field angular velocity and the magnetospheric 
electric current uj{^) and /(^) respectively These can 
only be determined when open field lines cross both sin- 
gular LS, the inner LS inside the ergosphere outside the 
black hole horizon, and the outer LS which is a deformed 
light cylinder. We developed an iterative numerical pro- 
cedure analogous to the one we implemented in the study 
of the pulsar magnetosphere. Iterating with respect to 
two functions simultaneously is a very difficult task and 
our method may certainly be improved in the future. 
Yet, we are confident enough to suggest that the dis- 
tributions are unique for the particular boundary con- 
ditions of the problem (all magnetic field lines thread 
the black hole horizon, all are open, and the surrounding 
disk is infinitely thin) since different initial conditions for 
5*, Lu and / evolve towards the same final solution (for 
each value of a that is). Our solutions may serve as test 
cases for time-dependent numerical codes. We plan to 
apply our method in the study of various astrophysically 
interesting situations where the surrounding disk is not 
infinitely thin but expands vertically with distance. 

All field lines that cross the black hole horizon are 
bound to cross the inner LS. This is not the case, though, 
for the outer LS. In fact, several numerical simulations 
exist in the literature where the black hole is threaded by 
a vertical uniform magnetic field, and yet they converge 
to certain distributions of w(5') and /(5') (e.g. Komis- 
sarov & McKinney 2007, Palenzuela et al. 2011). We 
expect that in that case, the solution is mathematically 
indeterminate, and the solutions shown depend on the 
particular choice of outer numerical boundary conditions. 
For example, the lateral boundary conditions in Palen- 
zuela et al. (2011) are not perfectly absorbing but reflect 
outgoing waves as is manifested in the oscillations seen 
in their figure 2. We plan to study such magnetic field 
configurations in a future work and check whether we can 
indeed freely specify one of the two free functions, u;(^') 
or /(^). We will thus check whether the results of previ- 
ous time-dependent numerical simulations are unique or 
not. 

We would like to caution the reader about certain nu- 
merical complications that arise in all numerical simu- 
lations of the spinning black hole magnetosphere, time 
dependent and steady state. Numerical simulations that 
do not adequately resolve both the inner and outer LS 
should not be trusted. We faced a similar problem in 
our study of the pulsar magnetosphere. In that case, 
there exists no inner LS, and in order to study the tran- 
sition through the outer LS (the light cylinder) we were 



forced to implement our inner boundary condition at a 
radial distance more than a hundred times larger than 
the actual neutron star radius (Spitkovsky 2006, McK- 
inney 2006, Contopoulos & Kalapotharakos 2010). This, 
however, posed no problem since the angular velocity of 
the field lines is set by the rotation of the pulsar. The 
black hole problem is more complicated because the black 
hole surface is not infinitely conducting, and the numer- 
ical code must adequately resolve both the inner bound- 
ary and the light cylinder. Notice that our numerical 
method works best as a/M approaches unity since in 
that case the inner LS is well separated from the event 
horizon, and the outer LS is at its closest approach to the 
origin of the numerical grid where its radial resolution is 
the highest. In that case, a numerical resolution of 256 
radial grid points corresponds to 34 grid cells between the 
horizon and the inner LS, and 94 grid cells between the 
two LS along the equator. When a/M > 0.9, our results 
are rather resolution independent. On the other hand, 
as a/M approaches zero, the ergosphere shrinks in width 
and the inner LS approaches so close to the event hori- 
zon that our numerical grid does not have enough radial 
resolution for the iterative procedure that yields w(^') 
to properly worlfl. At the same time, as a/M decreases 
so does f^BH, pushing the outer LS so far out that our 
radially expanding numerical grid has very low radial res- 
olution at the position of the outer LS for the iterative 
procedure to properly converge. Similarly, the resolu- 
tion of our radially expanding grid on the light cylinder 
becomes worse as — >■ 0, and this too complicates the 
convergence of our iterative method. Other numerical 
approaches certainly face similar complications at short 
and/or at large distances. Our method may be improved 
by a reformulation in radially expanding cylindrical coor- 
dinates, in Kerr-Schild coordinates, or in standard spher- 
ical coordinates with a numerical integration outer radial 
boundary. 

In summary, the force-free spinning black hole magne- 
tosphere is not very different from the force-free axisym- 
metric pulsar magnetosphere. In fact, beyond a distance 
on the order of a few times the black hole radius, the two 
problems are almost identical. An important corollary of 
this similarity is that a spinning black hole surrounded by 
a thin equatorial plasma disk does not generate a colli- 
mated outflow, only an uncollimated wind (as in Con- 
topoulos, Kazanas & Fendt 1999). Thus, the problem 
of the black hole jet formation remains open. A possi- 
ble resolution must involve a surrounding thick disk, or 
a surrounding magnetic disk wind that gradually coUi- 
mate the initially monopolar black hole wind into a col- 
limated jet outflow, as seen in the numerical simulations 
of Tchekhovskoy, Narayan & McKinney (2010). It is in- 
teresting to note here, again similarly to pulsars, that 
beyond the light cylinder the Lorentz factor V in the 
monopolar wind increases linearly with cylindrical dis- 



^ The solutions for a/M = 0.7 & 0.8 were obtained with a radial 
grid resolution of 1024. For even smaller values of a, our numer- 
ical integration fails because the inner LS is practically indistin- 
guishable from the event horizon. We can bypass this problem by 
choosing the inner boundary condition some distance outside the 
black hole horizon, set ui Ri O.Sf^BHi and iterate only with respect 
to the distribution of poloidal electric current /('I') that allows for 
a smooth crossing of the outer LS. 
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tance as 



r cx 



rsin6' rr^BHsin^ 



c/uj 



2c 



(20) 



(Contopoulos & Kazanas 2002), possibly up to the fast 
magnetosonic point of the outflow. Therefore, accord- 
ing to Eq. pO|. collimated field/flow lines that do not 
extend as far from the axis as uncollimated ones are ex- 
pected to reach lower Lorentz factors. Of course, the 
final Lorentz factor of such outflows will also depend on 
the mass loading of a given magnetic field lines, but such 
considerations are beyond the scope of the present work. 
We plan to continue our investigation of the interrelation 
between the black hole and disk outflows in the future. 
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